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Abstract: The cultivated area in artificial oases is deeply influenced by global climate change and human 
activities. Thus, forecasting cultivated area in artificial oases under climate change and human activities is of 
great significance. In this study, an approach named GD-HM-PSWROAM, consisting of general circulation 
model downscaling (GD), hydrological model (HM), and planting structure and water resource optimal 
allocation model (PSWROAM), was developed and applied in the irrigation district of the Manas River 
Basin in Xinjiang Uygur Autonomous Region of China to forecast the cultivated area tendency. 
Furthermore, the catchment export of the MIKE11 HD/NAM model was set to the Kensiwate 
hydrological station. The results show that the downscaling effects of temperature can be fairly satisfying, 
while those of precipitation may be not satisfying but acceptable. Simulation capacity of the MIKE11 
HD/NAM model on the discharge in the Kensiwate hydrological station can meet the requirements of 
running the PSWROAM. The accuracy of the PSWROAM indicated that this model can perform well in 
predicting the change of cultivated area at the decadal scale. The cultivated area in the Manas River Basin 
under current human activities may be generally decreasing due to the climate change. But the adverse 
effects of climate change can be weakened or even eliminated through positive human activities. The 
cultivated area in the Manas River Basin may even be increasing under assumed human activities and future 
climate scenarios. The effects of human activities in the future can be generally predicted and quantified 
according to the cultivated area trends under current human activities and the situations in the study area. 
Overall, it is rational and acceptable to forecast the cultivated area tendency in artificial oases under future 
climate change and human activities through the GD-HM-PSWROAM approach. 
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1 Introduction 


In recent years, global warming has become a world-wide environmental issue. The global 
averaged temperature has increased by 0.85°C during the period 1880—2012 (IPCC, 2013). The 
average precipitation over the mid-latitude land areas of the Northern Hemisphere has increased 
since 1901 (IPCC, 2013). In arid and semi-arid regions, runoff which generally originates in 
mountainous regions and dissipate in piedmont oases is the lifeline for the artificial oases (Zhou, 
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1999; Foda et al., 2017; Oduor and Mabanga, 2018). Further, the population in catchments is 
usually sparse, and the influence of human activities on runoff generation is fairly little (Li et al., 
2007; Chen et al., 2013; Xu et al., 2013; Chen et al., 2015). Therefore, runoff yield and 
concentration in arid and semi-arid regions mainly depend on the natural conditions, such as 
climate and the types and amounts of water resource in catchments (Zhou, 1999). The main water 
resource types in most catchments in arid and semi-arid regions are glaciers, snow and precipitation 
(Allan, 1996; Zhou, 1999; Dawoud, 2013). Global climate change have marked impacts on glaciers, 
snow and runoff in catchments (Chen et al., 2007; Pelto, 2008; Wang et al., 2017). Many previous 
researches showed that the glaciers are melting and the equilibrium-line altitudes are increasing in 
many regions since the 20" century (Fujita, 2008; Racoviteanu et al., 2008; Bolch et al., 2010a, b; 
Bhambri et al., 2011; Lambrecht et al., 2011; Scherler et al., 2011; Beedle et al., 2014; Zemp et al., 
2015; Cook et al., 2016). Furthermore, the water equivalent decreased in some arid and semi-arid 
regions (Aizen et al., 1997). Changes of glaciers, snow and precipitation have led to the variations 
in runoff (Chiew and Mcmahon, 2002). Therefore, it is necessary to analyze the trends of runoff 
under climate change in the future to better manage water resource in catchments. 

The artificial oases are the gathering places for human activities. Taking up only 3%-5% of the 
land area in arid and semi-arid regions, the artificial oases feed 90% of the population and create 
95% of the total output value (Lai, 2005). The main industry in artificial oases is agriculture. Water 
shortage is a very serious problem in most artificial oases. The percentage of water consumption 
for agriculture is 70% at the global scale, which is 85%—90% higher than that in water-deficient 
countries such as Asia and Africa (Delgado and Rodriguez, 2014). In the North Africa and Middle 
East, water shortage has resulted in health issues and even conflicts between countries (Shuval and 
Isaac, 1994; Aeschbacher et al., 2005; Johnson et al., 2016; Qian et al., 2016). In China, the water 
resource utilization efficiency is higher than 90% in most artificial oases, and 90% of water 
consumption is used in agriculture (Zhang et al., 2014). Many reservoirs in mountain passes are 
being constructed and canals and channels are being water-saving repaired to relieve the acute 
water shortage in artificial oases. Beyond that, water-saving irrigation technologies such as drip 
irrigation and sprinkling irrigation are widely used in many artificial oases in arid and semi-arid 
regions. 

Global warming has led to changes of temperature, precipitation and evaporation in arid and 
semi-arid regions, which will inevitably influence the river runoff in artificial oases and impact the 
available water resource and water utilization processes in the future. All of these changes are 
closely related to the cultivated area, which is of great significance to the economic growth, social 
security and ecological protection in artificial oases. However, few studies are focused on 
forecasting cultivated area in artificial oases under climate change and human activities. In this 
study, a cultivated area predicting approach was developed to explore the variations of cultivated 
area in artificial oases under future climate scenarios and human activities. The approach named 
GD-HM-PSWROAM for short is consisting of general circulation model downscaling (GD), 
hydrological model (HM), and planting structure and water resource optimal allocation model 
(PSWROAM). The methodology and performance of the approach was illustrated in details in the 
irrigation district of the Manas River Basin. 


2 Materials and methods 


2.1 Study area 


The irrigation district of the Manas River Basin is located in Xinjiang Uygur Autonomous Region, 
China. It is situated in the north of the Tianshan Mountains and south of the Gurbantunggut Desert 
(Fig. 1). The climate in the basin is of typical continental arid, with hot and dry summers, and 
chilly winters. The mean annual precipitation varies between 110 and 190 mm from the plains to 
the mountains (Ling et al., 2012). The annual mean temperature is about 6.6°C, and the mean 
annual evaporation ranges from 1500 to 2000 mm (Ling et al., 2012). The region is consisted of 
three sub-irrigation districts named the Shihezi, Xiayedi and Mosuowan. The mean cultivated area 
in the three sub-irrigation districts during the period 2001-2010 was 303.87, 626.47 and 506.20 
km, respectively (Shihezi Statistics Bureau, 2010). There are four plain reservoirs named the 
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Moguhu, Daquangou, Jiahezi and Yuejin in the irrigation district of the Manas River Basin. There 
is a mountain reservoir named the Kensiwate in the upstream of the Manas River. The effective 
capacity of the four plain reservoirs and Kensiwate reservoir is 1.73<10%, 0.36x10°, 0.65x108, 
0.75108 and 1.12x108 m°, respectively (Shihezi Statistics Bureau, 2015). 
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Fig. 1 Locations of the three sub-irrigation districts (Xiayedi, Mosuowan and Shihezi), river and canal systems, 
reservoirs (Moguhu, Daquangou, Jiahezi and Yuejin), Kensiwate hydrological station, Shihezi meteorology station 
and the selected European Centre for Medium-Range Weather Forecasts re-analysis Interim (ERA-Interim) points in 
the Manas River Basin 


The main water resources in the irrigation district of the Manas River Basin are river runoff and 
groundwater. The Manas River originates from the northern slope of the middle Tianshan 
Mountains. The annual runoff of the Manas River is 1.32*10° m°, with 69.6% of which appearing 
in June, July and August (Zhou, 1999). There are two hydrological stations named the Kensiwate 
hydrological station and Hongshanzui hydrological station in the main stream of the Manas River. 
The Hongshanzui hydrological station is located in the mountain pass, however, it was abandoned in 
1995. The Kensiwate hydrological station (43°58’00"N, 85°57'20"E; 900 m a.s.1.) is about 30 km 
upstream of the Hongshanzui hydrological station. About 95.5% of the runoff in the Manas River 
Basin can be monitored in the Kensiwate hydrological station. The annual exploitable groundwater 
of the Shihezi, Xiayedi and Mosuowan sub-irrigation districts is 1.96x10°, 0.78* 10% and 0.66x 10° 
m°, respectively (Shihezi Statistics Bureau, 2010). 


2.2 Data collection 


As the GD-HM-PSWROAM is an approach rather than a fixed model, some datasets for running it 
should be changed according to the selected models. Besides, these datasets have little effect on 
tendency of cultivated area. Therefore, the main datasets were induced in detail while other less 
important datasets were introduced in brief. 

Meteorological stations in catchments are fairly few in arid and semi-arid regions like Manas 
River Basin, due to high and steep terrain, and backward economic and technological situations. 
Furthermore, most of the stations are located in the mountain passes with discontinuous dataset. 
Even in the irrigation district, meteorological stations with observed temperature and precipitation 
for more than 30 years are few. Therefore, we selected the European Centre for Medium-Range 
Weather Forecasts re-analysis Interim (ERA-Interim) dataset to provide monthly mean daily 
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temperature (T), monthly mean daily maximum temperature (Tmax), monthly mean daily minimum 
temperature (Tmin) and monthly precipitation (P) as the substitution to the study due to its good 
performance (Hodges et al., 2010; Sylla et al., 2010; Dee et al., 2011; Weedon et al., 2014). Its 
spatial and temporal resolutions are 0.125°x0.125° and 1 month, respectively. The period of the 
dataset is from January 1979 to December 2012. The selected ERA-Interim points are shown in 
Figure 1. 

We used the general circulation models (GCMs), a kind of Navier-Stokes Equations basing on 
the rotating globe, to describe the physical processes of the atmosphere, ocean and land in the 
Earth (IPCC, 2013). The Representative Concentration Pathways (RCPs) are a set of scenarios 
containing the time series of all greenhouse gases, aero colloids and chemically active gases 
emissions and concentrations, and the scenarios of land use and cover (IPCC, 2013; Inatsu et al., 
2015). RCP2.6, RCP4.5 and RCP8.5 are the three commonly used RCPs in many GCMs, which 
represent the low, middle and high emission scenarios, respectively (IPCC, 2013; Inatsu et al., 
2015; Zhang et al., 2017). In this study, we obtained the outputs under RCP2.6, RCP4.5 and 
RCP8.5 scenarios of 28 GCMs (Table 1) by the Coupled Model Intercomparison Project, Phase 5 
from the node DEUTSCHES KLIMA RECHEM ZENTRUM (DKRZ,; https://www.dkrz.de/). 

The monthly mean runoff dataset during the period 1979-2012 in the Kensiwate hydrological 
station was provided by the Shihezi Water Conservancy Bureau. Other datasets such as crop 
irrigation schedule, water storage of reservoirs, and water surface evaporation used for the 
PSWROAM were obtained from the Statistical Yearbook of Shihezi (Shihezi Statistics Bureau, 
2010, 2015). 


Table 1 Spatial resolutions of the 28 GCMs provided by the DKRZ 


Model Spatial resolution (longitudexlatitude) Model Spatial resolution (longitude~latitude) 
BCC-CSMI-1 2.813°x2.791° GISS-E2-H 2.500°x2.000° 
BCC-CSM1-1-m 1.125°x1.121° GISS-E2-R 2.500°x2.000° 
BNU-ESM 2.813°x2.791° HadGEM2-AO 1.875°x1.250° 
CanESM2 2.813°x2.791° HadGEM2-ES 1.875°x1.250° 
CCSM4 1.250°x0.942° IPSL-CMS5A-LR 3.750°x 1.895° 
CESM1-CAM5 1.250°x0.942° IPSL-CM5A-MR 2.500°x1.268° 
CESM1-WACCM 2.500°x1.895° MIROC-ESM-CHEM 2.813°x2.791° 
CNRM-CMS5 1.406°x1.401° MIROC-ESM 2.813°x2.791° 
CSIRO-Mk3-6-0 1.875°x1.865° MIROCS5 1.406°x1.401° 
FGOALS-g2 2.813°x2.791° MPI-ESM-LR 1.875°x1.865° 
FIO-ESM 2.813°x2.791° MPI-ESM-MR 1.875°x1.865° 
GFDL-CM3 2.500°x2.000° MRI-CGCM3 1.125°x1.121° 
GFDL-ESM2G 2.500°x2.022° NorESM1-M 2.500°1.895° 
GFDL-ESM2M 2.500°x2.022° NorESM1-ME 2.500°x 1.895° 


Note: GCMs, general circulation models; DKRZ, DEUTSCHES KLIMA RECHEM ZENTRUM. 


2.3 Climate scene forecasting through GCM downscaling 


GCMs have become the most reliable and applicable approach to forecast climate scene in current 
studies. However, their spatial resolutions are too coarse and the biases are non-negligible for 
regional studies (Hewitson and Crane, 2006; Maurer and Hidalgo, 2008; Eden and Widmann, 2014; 
Inatsu et al., 2015). The downscaling approach could transform the coarse-scale outputs of GCMs 
to regional scales and correct the biases, which plays an important role in regional researches (Vrac 
et al., 2007; Inatsu et al., 2015). In this study, a combined approach of canonical correlation 
analysis (CCA) filtering, multi-model ensemble and extreme learning machine regressions 
(abbreviated as CEE) was applied in GCM downscaling because of its good performance (Zhang et 
al., 2017). The 7, Tmax, Tmin and P were downscaled in this study. 


JOURNAL OF ARID LAND 


In CEE downscaling, the CCA filtering was applied firstly to obtain the GCMs with the best 
simulation capacities on the T, Tmax, Tmin and P in the Manas River Basin. The T, Tmax, Tmin and P 
of the selected ERA-Interim points were extracted to constitute the ERA-Interim matrixes. Those 
of the nearest 25 points to each ERA-Interim point were extracted from the certain GCM outputs. 
Then, the corresponding GCM matrixes were constituted after removing redundancies. Finally, the 
CCA filtering between ERA-Interim matrixes and GCM matrixes was carried out, and the GCMs 
with the best simulation capacities were obtained. In this study, the number of the selected GCMs 
for T, Tmin, Tmax or P was 4 and the period of the ERA-Interim matrixes and GCM matrixes was 
from January 1979 to December 2012. 

After the CCA filtering, the multi-model ensembles of T, Tmax, Tmin of P were constituted by the 
4 selected GCMs. The weights of the 4 GCMs were equivalent here. However, the weight of the 
time series in every point of every GCM needs to be adjusted in the extreme learning machine 
(ELM) regressions. 

Finally, the spatial downscaling was carried out through ELM regressions. In this study, the ELM 
regressions in T, Tmax, Tmin and P downscaling were implemented from one ERA-Interim point to 
another. The corresponding dependent variable matrix was the ERA-Interim time series of T, Tmax, 
Tmin and P in the point. The corresponding independent variable matrix was constituted by the time 
series in the nearest 25 points of every selected GCM to the ERA-Interim point. It means that more 
than 100 time series could provide information. The calibration and validation periods in the ELM 
regressions was from January 1979 to December 2009 and from January 2010 to December 2012, 
respectively. 

After the GCM downscaling, the 7, Tmax, Tmin and P with needed resolution in the Manas River 
Basin were obtained. The potential evaporation (E) can be calculated through Har Formula (Eq. 1). 
Though the Penman Formula is more reliable, simulation capacities of most GCMs on wind speed 
are not reliable, and simulation results on sunshine duration cannot be obtained. Moreover, 
calculating Æ with more factors will increase the uncertainties of the results. Therefore, Har 
Formula rather than Penman Formula was used in this study. Moreover, the other parameters 
including the relative distance between solar and Earth, declination, sunset hour angle, and solar 
radiation were calculated by Equations 2—5, respectively. 

1 


E =ax Rax (Tnax — Tmin)? x (T +17.8)+ b, (1) 
dr =14.0.033xcos| 2x 2x], (2) 
365 

J 
= 0400xsin 2x2 -1.39), (3) 

365 
qs = arccos(—tang x tand) , (4) 
Ra = 24x 60x G, x dr x (@s x sin g x sin ô + cos Q x cos 6 x sin @s)/ T, (5) 


where Æ is the potential evaporation (mm); a and b are constants for calibration; Ra is the solar 
radiation absorbed by the top of atmosphere (MJ/(m7-d)); Tmax is the monthly mean daily maximum 
temperature (°C); Tmin is the monthly mean daily minimum temperature (°C); T is the mean daily 
temperature (°C); dr is the relative distance between solar and Earth; J is the day series; ô is the 
declination (rad); ws is the sunset hour angle (rad); ø is the latitude (rad); and Gs. is the solar 
constant (=0.0820 MJ/(m?-min)). 

As it is impossible to monitor the potential evaporation, we selected the Global Land-surface 
Evaporation (GLEAM, the Amsterdam Methodology) dataset (Miralles et al., 2011la) as the 
substitute dataset in calibrating the constants a and b in Har Formula. GLEAM is based on the 
Priestley and Taylor evaporation formula and the Gash analytical model of forest rainfall 
interception (Gash, 1979; Valente et al., 1997). It has been successfully validated over different 
climate conditions (Miralles et al., 2011a, b, c; Poveda and Zapata, 2016). The GLEAM dataset 
was adjusted referring to the observed water surface evaporation in the Kensiwate hydrological 
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station and Shihezi meteorological station before it was applied in calibrating the constants. 
2.4 Runoff forecasting through hydrological model 


Hydrological models have been widely applied in hydrological system simulation (Abbott et al., 
1986; Xu et al., 1996; Santhi et al., 2001; Mcmichael et al., 2006; Sahoo et al., 2006; Arnold et al., 
2012). MIKE11, developed by Danish Hydraulic Institute (DHI), is a module for simulating the 
runoff of rivers and canals, water quality and sediment transport (Shamsudin and Hashim, 2002). 
MIKE11 NAM is a rainfall-runoff model of MIKE11 RR (rainfall-runoff) module. NAM, one of 
the lumped and conceptual hydrological models advanced by Nielsen and Hansen (Makungo et al., 
2010), is mainly applied in simulating the rainfall-runoff process in natural watershed. The 
temperature, precipitation, potential evaporation and initial snow water equivalent datasets are 
necessary in simulating the runoff yield and concentration process in NAM. MIKE11 HD/NAM 
has been widely applied in many researches due to its good performance and less required datasets 
(e.g., Delphi, 2011; Doulgeris et al., 2012; Luo et al., 2013; Lin et al., 2014; Vazifedoost et al., 
2014; Gao et al., 2017). Therefore, in this study, we selected this module for the runoff simulation 
of the Manas River. 

About 95.5% of the runoff in the Manas River can be monitored in the Kensiwate hydrological 
station. In addition, the Kensiwate reservoir locates in only 2 km upstream of the Kensiwate 
hydrological station, and there are no tributaries between them. Selecting the Kensiwate 
hydrological station as the export means that the predicted runoff can be applied in the PPWROAM 
without adjustment. Therefore, the catchment export of the MIKE11 HD/NAM model was set to 
the Kensiwate hydrological station. 

The catchment boundary of the Kensiwate hydrological station (Fig. 2) and rivers inside the 
catchment were extracted from the Advanced Spaceborne Thermal Emission and Reflection 
Radiometer Global Digital Elevation Model (ASTER GDEM DEM; http://www.gscloud.cn/). The 
temperature and precipitation time series of six sub-catchments not shown were calculated from the 
ERA-Interim dataset through the Thiessen polygon approach. The potential evaporation time series 
was calculated from the temperature series through Har Formula. The initial snow water equivalent 
was extracted from the ERA-Interim dataset and the initial glacier water equivalent was obtained 
from the China Glacier Information System (Wu and Li, 2003). We set the calibration period and 
validation period from January 1979 to December 2009 and from January 2010 to December 2012, 
respectively, based on the period (January 1979 to December 2012) of the obtained runoff dataset 
in the Kensiwate hydrological station. 


2.5 Cultivated area forecasting through PPWROAM 


As more than 90% of the water consumption is for agriculture in many artificial oases in global, the 
cultivated area mainly depends on the water supply and it can be forecasted through PSWROAM. 
The PSWROAM is consist of the planting structure optimization model and water resource optimal 
allocation model. In most artificial oases, water utilization is considered as a constraint in building 
the model. In fact, the goal of water resource optimal allocation model is to guarantee the water 
supply for the planting structure optimization model. In this study, we constructed the water 
resource optimal allocation model according to the river and canal systems, reservoirs and 
sub-irrigation districts (Fig. 3). 

The mathematical functions of PSWROAM were generalized as follows. The objective function 
is described as Equation 6. 


3 >on 
max(OV) = nanl S35, x OV, , ) ; (6) 
i=l j=l 

where OV is the output value of the planting industry in the irrigation district of the Manas River 
Basin (USD); Si; is the cultivated area of the crop type numbered j in the sub-irrigation district 
numbered i (m°); OV;; is the output value per unit area of the crop type numbered j in the 
sub-irrigation district numbered i (USD); and n is the number of crop types in the irrigation district 
of the Manas River Basin. 
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The constraint of cultivated area for crops can be expressed as: 
S S072 8 o (7) 


i,j,min 
where Sjjmin is the minimum cultivated area of the crop type numbered j in the sub-irrigation 
district numbered i (m°); and S;;max is the maximum cultivated area of the crop type numbered j in 
the sub-irrigation district numbered i (m7). 

Moreover, the constraint of cultivated area for the sub-irrigation districts can be expressed as 
follows: 


S; min < S; < S; max 5) (8) 
where Simin is the minimum cultivated area of the sub-irrigation district numbered i (m°); S; is the 
cultivated area of the sub-irrigation district numbered i (m°); and Simax is the maximum cultivated 
area of the sub-irrigation district numbered i (mô). 


The constraint of output value for the sub-irrigation districts is shown in Equation 9. 


S OV, = YG, x OV, ,) > (9) 


j=l 


OV, min 
where OVimin is the minimum output value of the sub-irrigation district numbered i (USD); and OV; 
is the output value of the sub-irrigation district numbered 7 (USD). 

The constraint of water supply for the sub-irrigation districts can be expressed as: 


n 
Vki + Vione s DA x 1; + V k gw + Vis 2 (10) 
j=l 


where Vikir is the needed water quantity for agricultural irrigation of the sub-irrigation district 
numbered i in the month numbered k (m°); Vikotne is the needed water quantity for other industries 
of the sub-irrigation district numbered i in the month numbered k (m°); V;+ is the diverted water 
quantity of the canal numbered j flowing into the sub-irrigation district in the month numbered k 
(m°); n; is the conveyance efficiency of the canal numbered j; Vikgw is the groundwater exploitation 
of the sub-irrigation district numbered i in the month numbered k (m°); and Vikotav is the other 
available water quantity of the sub-irrigation district numbered i in the month numbered k (m°). 
The constraints of groundwater exploitation are shown in Equations 11 and 12: 
y, <y, <V, 


i,k,gw,min i,k,gw i,k,gw,max ? (l 1) 


12 


4 


aw S Vigwma > (12) 
k=1 
where Vik gw,min is the minimum groundwater exploitation of the sub-irrigation district numbered i 
in the month numbered k (m°); Vikgwmax is the maximum groundwater exploitation of the 
sub-irrigation district numbered i in the month numbered k (m°); and Vigwmax is the maximum 
annual groundwater exploitation of the sub-irrigation district numbered i (m°). 


The constraint of ecological water supply can be expressed as follows: 
Vrne S Vase X Moa + Vas X The» (13) 


k,ne 


where Vine is the needed ecological water quantity of the Manas River in the downstream of the 
Jiahezi reservoir in the month numbered k (m°); V24, and V6, are the diverted water quantity of the 
canal numbered 24 and 26 flowing into the sub-irrigation district in the month numbered k (m°), 
respectively; and 424 and 7426 are the conveyance efficiencies of the canal numbered 24 and 26, 
respectively. 

Water balance constraint of the points where canals join or separate is shown in Equation 14. 


Sg = aks (14) 
j=l j=l 


where V;kin is the diverted water quantity of the canal numbered j flowing into the point in the 
month numbered k (m°); and V; xt is the diverted water quantity of the canal numbered j flowing 
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from the point in the month numbered k (m°). 
Water balance constraint of reservoirs is expressed as: 


Via 


+y, 


j,k,ca,in 


y 


i,k,ev 


y 


i,k,le 


x7, 4 y, 


i,k,ot,in 


=V (15) 
j=l 

where Vix is the water storage of the reservoir numbered i in the end of the month numbered k-1 
(m°); Vj4ca,in is the diverted water quantity of the canal numbered j flowing into the reservoir in the 
month numbered k (m°); Vikotin is the other water quantity flowing into the reservoir numbered i in 
the month numbered k (m°); Vikev is the evaporation loss of the reservoir numbered i in the month 
numbered k (m°); Vikile is the leakage loss of the reservoir numbered i in the month numbered k 
(m°); V; kcat is the diverted water quantity of the canal numbered j' flowing from the reservoir in 
the month numbered k (m°); and V;x is the water storage of the reservoir numbered i in the end of 
the month numbered k (m°). 

Water storage constraint of reservoirs is expressed as Equation 16. 


V pais < i,k < V iieii , (16) 


where Vikmin is the minimum water storage of the reservoir numbered i in the end of the month 
numbered k (m°); and Vikmax is the maximum water storage of the reservoir numbered i in the end 
of the month numbered k (m°). 

The constraint of diverted water quantity for canals is shown in Equation 17. 


0<V..<V,., (17) 
j,k j,k,max 


where V; is the diverted water quantity of the canal numbered j in the month numbered k (m°) and 
V;kmax is the maximum diverted water quantity of the canal numbered j in the month numbered k 
(m°). 
Nonnegative constraint is shown as follows: 
0<x, (18) 
where x refers all the physical quantities in the model. 
Evaporation loss of reservoirs can be described as Equation 19: 


Viney = Eig X Sips (19) 


i,k,ev 
where Virev is the evaporation loss of the reservoir numbered i in the month numbered k (m°); Eix 
is the water surface evaporation per unit area of the reservoir numbered i in the month numbered k 
(m?/m’); and Six is the average water surface area of the reservoir numbered i in the month 
numbered k (m°). 
Moreover, leakage loss of reservoirs is shown in Equation 20: 


Vi ate =Q pX Vins + Vix) /2, (20) 


i 


where Vike is the leakage loss of the reservoir numbered i in the month numbered k (m°); aix is the 
leakage coefficient of the reservoir numbered i in the month numbered k (m°); Vix- is the water 
storage of the reservoir numbered i in the end of the month numbered k—1 (m°); and V;,; is the water 
storage of the reservoir numbered i in the end of the month numbered k (m°). 


2.6 Evaluation indices 


The coefficient of determination (R°) and relative mean absolute residual (RMAR) were selected to 
evaluate the performance of the ELM regressions (Zhang et al., 2017), which can be calculated 
through Equations 21 and 22, respectively. The R?, RMAR and Nash coefficient (Rys) were used to 
evaluate the performance of MIKE11 HD/NAM model. The Rys is expressed as Equation 23. 


R1- DOs at , (21) 
Yn = Vor) 
RMAR = > Yoo — Ysi (22) 


nx(max(y,,)—min(y,,))” 
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yg =1- See) (23) 
= S Oa = Va) Í 


where n is the length of the time series; yo is the observed time series; ysi is the simulated time 
series; and y,, is the mean of the observed time series. 


3 Results 


3.1 Performance of the GCM downscaling 


The mean R? and RMAR values of the GCM downscaling in the selected ERA-Interim points are 
listed in Table 2. The mean R? values of the ELM regressions on T, Tmax and Tmin were all >0.967, 
which means that the effect of CEE downscaling on temperature is satisfying. The mean RMAR 
values of the ELM regressions on T, Tmax and Tmin were all >0.038, indicating that the regression 
error is about 4.0% of the series variation. Compared with temperature, the performance of CEE 
downscaling on P is not satisfying because the mean R? values were lower (<0.405) and the mean 
RMAR values were higher (0.169-0.174) under three climate scenarios (RCP2.6, RCP4.5 and 
RCP8.5). However, the effect of CEE downscaling on precipitation is acceptable in consideration 
of the vast spatial and temporal variations of precipitation in the study area (Zhang et al., 2017). In 
the Manas River Basin, the snow/glacier meltwater takes up more than 70% of the runoff (Zhou, 
1999; Chen et al., 2015). Therefore, the influence of the downscaling effect of precipitation on 
runoff and cultivated area is small. 


Table 2 Mean R? and RMAR values of the ELM regressions in all selected ERA-Interim points 
Climate factor Scenario R? RMAR Climate factor Scenario R? RMAR 

RCP2.6 0.967 0.040 RCP2.6 0.967 0.043 
T RCP4.5 0.971 0.040 Tmin RCP4.5 0.971 0.041 
RCP8.5 0.973 0.038 RCP8.5 0.973 0.038 
RCP2.6 0.973 0.038 RCP2.6 0.378 0.172 
Tras RCP4.5 0.967 0.041 P RCP4.5 0.403 0.169 
RCP8.5 0.970 0.039 RCP8.5 0.405 0.174 


Note: R?, coefficient of determination; RMAR, relative mean absolute residual; ELM, extreme learning machine; T, mean daily 
temperature; Tmax, monthly mean daily maximum temperature; Tmin, monthly mean daily minimum temperature; P, monthly precipitation; 
RCP, Representative Concentration Pathway. 


3.2 Simulation capacity of the MIKE11 HD/NAM model 


Although the calibration period of the MIKE11 HD/NAM model was from January 1979 to 
December 2009, evaluating indices were calculated during the period from January 1980 to 
December 2009 because the performance of the model in 1979 was greatly affected by the initial 
conditions such as initial discharge or ground water level (Delphi, 2011; Doulgeris et al., 2012; 
Vazifedoost et al., 2014). Indices for the validation were calculated in the period from January 2010 
to December 2012. The indices used for evaluating the MIKE11 HD/NAM model are listed in 
Table 3, and the simulated and observed runoff in the Kensiwate hydrological station is shown in 
Figure 4. As shown in Table 3, the R? and Rys values during the calibration and validation periods 
were >0.893. It means that the MIKE11 HD/NAM model is performed well on the simulation of 
the runoff in the Kensiwate hydrological station. Furthermore, the RMAR during the validation 
period was 0.056, showing that the simulation error is about 5.6% for the maximum runoff 
variation in the Kensiwate hydrological station. Therefore, the MIKE11 HD/NAM model is 
satisfactory for running the PSWROAM. 


3.3 Simulation capacity of PSWROAM in predicting cultivated area 


It is difficult to verify the accuracy of the PSWROAM from 1979 year by year since many datasets 
are required to drive the model. Moreover, the effective capacities of reservoirs are always changing 
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Table 3 Indices used for evaluating the MIKE11 HD/NAM model during the calibration and validation periods 
Period R Rys RMAR 
Calibration period 0.895 0.893 0.033 
Validation period 0.921 0.921 0.056 
Note: Rys, Nash coefficient. 
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Fig. 4 Simulated and observed runoff in January in the Kensiwate hydrological station during the period 1979— 
2012 


rapidly because the sedimentation and the corresponding monitoring dataset are confidential. 
Therefore, the PSWROAM was only used to simulate the cultivated area at the decadal scale. It 
was verified during the periods 1991-1995, 1996-2000, 2001-2005 and 2006-2010, as shown in 
Table 4. The simulated cultivated area was fairly close to the actual cultivated area, with the R? 
value of 0.97. Therefore, the accuracy of the PSWROAM is satisfactory and it can be used to 
analyze the variations of cultivated area under the future climate scenarios at the decadal scale. 


Table 4 Simulated and actual cultivated area in the irrigation district of the Manas River Basin 


Cultivated area 1991-1995 1996-2000 2001-2005 2006-2010 
Simulated cultivated area (km?) 1281.47 1336.74 1435.06 1436.53 
Actual cultivated area (km?) 1298.56 1320.48 1430.23 1442.84 


3.4 Trends of temperature, precipitation and potential evaporation under climate scenarios 


After the CEE downscaling, the monthly mean potential evaporation was calculated by Har 
Formula. Then, the decadal annual mean temperature and mean annual precipitation and potential 
evaporation in the catchment of the Kensiwate hydrological station and in the irrigation district of 
the Manas River Basin were calculated. The results are shown in Figures 5-7, respectively. 
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Fig. 5 Decadal annual mean temperature in the catchment of the Kensiwate hydrological station (C) and in the 
irrigation district (I) of the Manas River Basin under RCP2.6, RCP4.5 and RCP8.5 scenarios during the period 
2001-2100. Values in the brackets are the variation rates at the decadal scale. 
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Fig. 6 Decadal mean annual precipitation in the catchment of the Kensiwate hydrological station (C) and in the 
irrigation district (I) of the Manas River Basin under RCP2.6, RCP4.5 and RCP8.5 scenarios during the period 
2001-2100. Values in the brackets are the variation rates at the decadal scale. 
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Fig. 7 Decadal mean annual potential evaporation in the catchment of the Kensiwate hydrological station (C) and 
in the irrigation district (I) of the Manas River Basin under RCP2.6, RCP4.5 and RCP8.5 scenarios during the 
period 2001-2100. Values in the brackets are the variation rates at the decadal scale. 


The decadal annual mean temperature in the catchment will increase under RCP2.6, RCP4.5 and 
RCP8.5 scenario (Fig. 5), which may speed up the glacier melting. As shown in Figure 6, the 
decadal mean annual precipitation in the catchment will generally increase. Compared with the 
period 2001-2010, the decadal mean annual potential evaporation of the catchment during the 
period 2091-2100 will increase under RCP2.6 and RCP4.5 scenarios but decrease under RCP8.5 
scenario (Fig. 7). 

The decadal annual mean temperature in the irrigation district of the Manas River Basin will 
increase under RCP2.6, RCP4.5 and RCP8.5 scenarios, with the increasing rate being much higher 
than that in the catchment of the Kensiwate hydrological station (Fig. 5). The decadal mean annual 
precipitation will increase under RCP2.6 scenario but decrease under RCP4.5 and RCP8.5 
scenarios in the irrigation district (Fig. 6). Furthermore, the decadal mean annual potential 
evaporation in the irrigation district of the Manas River Basin will increase under RCP2.6 and 
RCP4.5 scenarios but decrease under RCP8.5 scenario (Fig. 7). 


3.5 Trends of annual runoff in the Kensiwate hydrological station under climate scenarios 


The temperature, precipitation and potential evaporation of sub-catchments during the period 
2013-2100 were calculated through the Thiessen Polygon approach. After that, the MIKE11 
HD/NAM model was driven and the annual runoff and corresponding glacier and snow water 
equivalent during this period were obtained. Since the inter-annual variations of annual runoff were 
intense and the PSWROAM was driven at the decadal scale, the decadal mean annual runoff was 
obtained during the period 2013—2100 under the three scenarios and the results are shown in Figure 
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8. Moreover, the trends of mean glacier and snow water equivalent in the catchment under the three 
scenarios were calculated (Fig. 9) because this water equivalent is closely related to the annual 
runoff in the Kensiwate hydrological station. 

As shown in Figure 8, the annual runoff in the Kensiwate hydrological station will generally 
decrease under RCP2.6 and RCP4.5 scenarios but increase under RCP8.5 scenario. The decadal 
mean annual runoff during the period 2091-2100 will be 3.70x10’-5.50x10’ m? higher than that 
during the period 2001-2010. As shown in Figure 9, the glacier and snow water equivalent in the 
catchment of the Kensiwate hydrological station will decrease rapidly under the three climate 
scenarios. According to the results of MIKE11 HD/NAM model in this study and previous study 
(Zhou, 1999), glacier and snow meltwater contributes about 70% of the total runoff. Therefore, the 
increase of runoff is mainly attributed to the melting of glacier and snow. The basic reason for the 
decrease of runoff under RCP2.6, RCP4.5 and RCP8.5 scenarios is the decrease of the glacier and 
snow water equivalent. 
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Fig. 8 Decadal mean annual runoff in the Kensiwate hydrological station under RCP2.6, RCP4.5 and RCP8.5 
scenarios during the period 2001—2100. Values in the brackets are the variation rates at the decadal scale. 
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Fig. 9 Decadal mean glacier and snow water equivalent in the catchment of the Kensiwate hydrological station 
under RCP2.6, RCP4.5 and RCP8.5 scenarios during the period 2001-2100. Values in the brackets are the variation 
rates at the decadal scale. 


3.6 Trends of cultivated area under climate scenarios and current human activities 


Due to water shortage in the Xiayedi and Mosuowan sub-irrigation districts is serious and there is 
no more land to expand farmland in the Shihezi sub-irrigation district, the maximum cultivated area 
in the irrigation district of the Manas River Basin under current human activities was set to equal to 
that during the period 2001-2010. Finally, the PSWROAM was driven and the cultivated area in 
the irrigation district of the Manas River Basin under the three climate scenarios and current human 
activities was predicted (Fig. 10). 
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The cultivated area in the irrigation district of the Manas River Basin under the three climate 
scenarios and current human activities will generally decrease. The decreasing rate of cultivated 
area under RCP2.6, RCP4.5 and RCP8.5 is 15.796, 8.193 and 2.990 km?/10a, respectively. The 
decadal mean cultivated area during the period 2091-2100 will be 0.000—122.900 km? lower than 
that during the period 2001-2010. 
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Fig. 10 Decadal mean cultivated area in the irrigation district of the Manas River Basin under current human 
activities and RCP2.6, RCP4.5 and RCP8.5 scenarios during the period 2001—2100. Values in the brackets are the 
variation rates at the decadal scale. 


4 Discussion 


4.1 Possible human activities to cope with the cultivated area decrease and the trends of 
cultivated area trends under assumed human activities 


As above-mentioned, cultivated area in the irrigation district of the Manas River Basin will 
decrease under current human activities and climate change. However, the decrease in cultivated 
area indicates land degradation in this region, which is very dangerous for the ecological 
environment and economic development. Therefore, measures should be taken to cope with this 
situation. The possible measures may be included: improving the conveyance efficiency of the 
canals, reducing the crop water consumption and abandoning the plain reservoirs. 

The conveyance efficiency of the canal system in the Shihezi, Xiayedi and Mosuowan 
sub-irrigation districts is 0.82, 0.77 and 0.77, respectively (Shihezi Statistics Bureau, 2010), which 
can be improved in the future. The conveyance efficiency of the canal system in the three 
sub-irrigation districts was assumed to increase at the rate of 0.025/10a since 2021, until 0.900. As 
the conveyance efficiency is lower than 1.000, the assumption is feasible. 

Even though the water requirement of cropping systems in the irrigation district of the Manas 
River Basin is relatively low, it is possible to reduce the water requirement by considering the 
technology progress such as deficit irrigation (Romero, et al., 2004). Therefore, the monthly water 
requirement of cropping systems in the three sub-irrigation districts was assumed to decrease at the 
rate of 0.15 m?/(hm?-10a) on the base of that under current human activities. 

After the operation of the Kensiwate reservoir in 2015, the flood prevention function of the 
Jiahezi reservoir has been taken over by the Kensiwate reservoir. The sedimentation capacity of the 
Jiahezi reservoir has reached 6.6010’ m?, which is very serious because its total capacity is only 
11.00x10’ m°. Besides, the Jiahezi reservoir is a defective reservoir. Therefore, the Jiahezi reservoir 
may be abandoned in the future. In this paper, the Jiahezi reservoir was assumed to be abandoned 
in 2031 (Yang et al., 2013). 

The water resource may be sufficient if the measures listed above are taken. Due to the economic 
growth and human survival stress, the cultivated area may increase if water resource is sufficient. 
However, there is no redundant land for farmland area in the Shihezi sub-irrigation district. The 
cultivated area expansion is only possible in the Xiayedi and Mosuowan sub-irrigation districts. 
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Therefore, the maximum cultivated area in the Xiayedi and Mosuowan sub-irrigation districts after 
2021 was assumed to be 1.5 times of that during the period 2001-2010. 

We adjusted parameters according to the quantified human activities assumed above, and drove 
the PSWROAM subsequently. Then the decadal mean cultivated area in the irrigation area of the 
Manas River Basin under assumed human activities was obtained (Fig. 11). 

The cultivated area of the Manas River Basin under assumed human activities will generally 
increase. It means that the adverse effect of climate change on cultivated area may be weakened or 
eliminated through positive human activities. Comparing to that under current human activities, the 
cultivated area under RCP2.6 and RCP8.5 scenarios is much higher during the period 2031-2040, 
which is mainly due to the abandonment of the Jiahezi reservoir. As shown in Table 5, abandoning 
the Jiahezi reservoir will obviously reduce the evaporation and leakage losses of the reservoir 
groups, which is fairly beneficial for the planting industry. However, the cultivated area increase 
after 2041 is mainly due to the improved conveyance efficiency of canals and reduced water 
requirement of crops. Overall, the decadal mean cultivated area during the period 2091-2100 will 
be 78.360—185.070 km? higher than that during 2001-2010. 
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Fig. 11 Decadal mean cultivated area of the Manas River Basin under assumed human activities and RCP2.6, 
RCP4.5 and RCP8.5 scenarios during the period 2001-2100. Values in the brackets are the variation rates at the 
decadal scale. 


Table 5 Evaporation and leakage losses under current and assumed human activities during the period 2031-2040 


; Evaporation loss (x 10’ m°) Leakage loss (x 10’ m°) 
Scenario 
Current Assumed Current Assumed 
RCP2.6 5.668 3.944 2.905 2.170 
RCP4.5 5.976 4.006 3.207 2.221 
RCP8.5 5.708 4.041 2.862 2.191 


4.2 Evaluation and application of GD-HM-PSWROAM approach 


The GD-HM-PSWROAM is consist of GCM downscaling, hydrological model, planting structure 
and water resource optimal allocation model. It provides an approach rather than a fixed model. 
Therefore, its advantages and disadvantages were analyzed in view of an approach rather than a 
fixed model. 

The GD-HM-PSWROAM provides an approach in forecasting the cultivated area in artificial 
oases under climate change and human activities. Its applicability in artificial oases was 
demonstrated rational. The effects of GCM downscaling on temperature can be fairly satisfying, 
while its effect on precipitation is not satisfying but acceptable at least. Simulation capacity of the 
hydrological model can be satisfying for driving the PPWROAM. The accuracy of the PPWROAM 
can be fairly satisfying for analyzing the cultivated area at the decadal scale. In general, accuracy of 
the GD-HM-PSWROAM in predicting cultivated area at the decadal scale is fairly satisfying. 

Unlike other fixed models, the specific models in GD-HM-PSWROAM need to be selected and 
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adjusted on the basis of the practical situations of the study regions. Its performance is 
tremendously influenced by errors and uncertainties of the input datasets and the selected specific 
models. Raising the required datasets of the GD-HM-PSWROAM is fairly intractable. The 
matching among the input datasets and outputs datasets of the selected models is fairly uneasy. 

The GD-HM-PSWROAM was only applied in artificial oases in this study. However, it does not 
mean that the approach cannot be applied in other regions. But there are a couple of points to be 
noticed when applying this approach in other regions. First, the monitoring datasets rather than the 
reanalysis datasets should be selected if the number and distribution of meteorological stations are 
satisfying. The monitoring datasets are more exact and reliable than the reanalysis datasets in 
general (Abbott et al., 1986; Hodges et al., 2010; Miralles et al., 2011c). Second, a distributed 
hydrological model rather than a lumped or semi-distributed model should be employed if the 
required datasets are satisfying. In lumped model and or semi-distributed model, the glacier and 
snow are usually assumed to be distributed evenly in the basins or sub-basins (Xu, 2009). However, 
the glacier and snow water equivalent can be precisely located into each cell in distributed models 
(Xu et al., 1996; Shamsudin and Hashim, 2002; Delphi, 2011; Luo et al., 2013; Lin et al., 2014; 
Vazifedoost et al., 2014). Third, a relatively closed region rather than an open region is 
recommended. The pass-by water resource in study regions will not predicted in the 
GD-HM-PSWROAM, which may lead to errors and uncertainties. In some open regions, the 
quantity of the pass-by water resource is even higher than that of the indigenous water resource. 
Finally, precisely forecasting and quantifying the human activities are of great importance. The 
regulation effects of hydrological structures such as the reservoirs, sluice gates and pumps are 
enormous. Constructing or abandoning the reservoirs and channels will alter the PPWROAM on 
some degree. Drawing water for other industries will slow down the water supply for the planting 
industries. Applying water-saving irrigation technologies may slash the water requirement of the 
planting industries. To sum up, the GD-HM-PSWROAM provides an approach rather than a fixed 
model. Selecting proper models in GD-HM-PSWROAM will lead to better results. In other words, 
the improper models may lead to errors and uncertainties. 


5 Conclusions 


In this paper, the GD-HM-PSWROAM approach was developed to forecast the cultivated area in 
artificial oases under climate change and human activities at the decadal scale. The 
GD-HM-PSWROAM is consisting of GCM downscaling, hydrological model, planting structure 
and water resource optimal allocation model. The performance and methodology of the 
GD-HM-PSWROAM was illustrated in details in the Manas River irrigation district Basin. The 
effect of CEE downscaling effect on temperature can be satisfying since the mean R? is 0.97 and 
the mean RMAR is 0.04. Although the effect of CEE downscaling on precipitation may be not 
satisfying but it is acceptable at least since the R? is 0.40 and the RMAR is 0.17. The simulation 
capacity of MIKE11 HD/NAM model on the discharge in the Kensiwate hydrological station can 
be satisfying for driving the PSWROAM since the R? and Rys values during the validation period 
are all 0.92. The accuracy of the PSWROAM is satisfying for analyzing the cultivated area at the 
decadal scale with the R? of 0.97. 

The cultivated area of the Manas River Basin under current human activities will generally 
decrease under RCP2.6, RCP4.5 and RCP8.5 scenarios. Measures need to be taken to cope with the 
decreasing trend. The adverse effects of climate change on cultivated area may be weakened or 
eliminated through positive human activities. 

Overall, we can predict the effects of human activities according to the trends of cultivated area 
under current human activities and situations in the study area. Then the trends of cultivated area 
under climate change and human activities can be forecasted through the PSWROAM. Therefore, it 
is rational to forecast the cultivated area in artificial oases under climate change and human 
activities through the GD-HM-PSWROAM approach. 
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